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Abstract. 

We present a theoretical analysis of the environment effects on charge transport in 
double-stranded synthetic poly(G)-poly(C) DNA molecules attached to two ideal leads. 
Coupling of the DNA to the environment results in two effects: (i) localization of carrier 
functions due to the static disorder and (ii) phonon-induced scattering of the carrier 
between these localized states, resulting in hopping conductivity. A nonlinear Pauli 
master equation for populations of localized states is used to describe the hopping 
transport and calculate the electric current as a function of the applied bias. We 
demonstrate that, although the electronic gap in the density of states shrinks as the 
disorder increases, the voltage gap in the I —V characteristics becomes wider. Simple 
physical explanation of this effect is provided. 
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1. Introduction 

Electronic transport tlirougli DNA molecules attached to leads still remains a 
controversial topic. A number of experiments on electrical transport through dry and 
wet DNA molecules revealed a variety of results. Double stranded DNA demonstrated 
proximity-induced superconducting [1], metallic [2-4], semiconducting [5-10] and 
insulating [11, 12] behaviors. The observed differences are often attributed to contact 
effects, coupling of the DNA to the environment, and to the sequence of nucleotides. 
Due to the diversity of experimental results, no consensus on mechanisms responsible 
for charge transport in the DNA has been achieved so far. 

Semiconducting behavior of double stranded synthetic poly(G)-poly(C) DNA was 
established experimentally. The I — V characteristic of the molecule revealed a voltage 
gap [5]. Effective Hamiltonian models, based on the tight-binding approximation [13- 
23], provided a reasonable description of the semiconductor gap observed in experiments 
with a minimum set of adjustable parameters. 

Less effort has been devoted to study the effects of molecular vibrations on the 
electric current through DNA. It was found that hopping of the charge between 
the sites of guanine (G) traps and the charge-phonon coupling results in a staircase 
structure of the I — V characteristics [24]. The influence of vibrational modes on the 
electronic properties of various types of DNA molecules (synthetic and natural) was 
addressed in Ref. [25]. It was argued that charge transport is dominated by quasi- 
ballistic contributions in homogeneous DNA and the zero-bias conductance is enhanced 
by the coupling to vibrations. Dissipative effects in the electronic transport through 
DNA molecular wires, comprising counter-ions and hydration shells, were investigated 
in Ref. [26]. A bath-induced pseudo-gap opens in the strong-coupling regime and a 
crossover from tunneling to phonon-assisted transport was observed with increasing 
temperature. Also, it was claimed that disorder effects smear the electronic band but 
have negligible impact on the formation of the pseudo-gap. 

In this work we focus on the effect of the environment on the charge transport 
properties of DNA molecules. To this end, we use the tight-binding ladder model [13] 
with random base-pair energies to describe the DNA electronic states, when disorder 
localizes electronic states. Charge transport is then mediated by phonon-assisted 
hopping between these states, which is described by means of a nonlinear Pauli master 
equation enabling us to calculate the current as a function of the applied voltage. We 
study the dependence of the I — V characteristics on the magnitude of disorder. 

2. Ladder model of DNA 

We consider the ladder-like model [13] of the poly(G)-poly(c) DNA based on a tight- 
binding Hamiltonian in the nearest-neighbor approximation. Figure [1] represents the 
schematics of the model. 
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Figure 1. Schematic view of a fragment of poly(G)-poly(C) DNA molecules, excluding 
the sugar-phosphate backbone, coupled to ideal leads. 

The corresponding Hamiltonian reads: 

- ^n) = tec {€+1 + + tcG€ , 

{E - 5^) = tcG {€+1 + €-l) + tcG€ ■ (1) 

Here superscripts G and C label a strand, and are energies of base molecules, 
tcG and tec are intra-strand transfer interactions (from now on we use tec = ice for 
simplicity), while tec is the inter-strand interaction. In the absence of disorder, site 
energies are the same along each strand. 

3. Coupling to the environment 

Interactions with a random environment of solute molecules and ions surrounding the 
DNA molecule can (i) result in perturbation of the energies of base molecules and (ii) give 
rise to coupling to phonons of the bath. We account for the former effect by considering 
random base energies uniformly distributed within a box of width A and centered around 
the unperturbed energies Sg or Sc, depending on the strand. Disorder, originated from 
e.g. electrostatic interactions with solute ions, can be large and all electronic states are 
therefore localized at segments which are typically shorter than the length of the DNA 
molecule. The system would therefore have exponentially small transmission coefficient 
and direct tunneling through the system is expected to be suppressed. On the other 
hand, electron-phonon coupling can result in phonon-assisted hopping between these 
localized states, i.e. incoherent charge transport. 

For each realization of disorder, we diagonalize the Hamiltonian ([1]) and calculate 
the scattering rate between the eigenstate ipp (with energy Ep) and another one ipa 
(with energy Ea) according to (see Refs. [27-31] for further details): 

W^p = S{Ep - E^) I^p E{Ep - E^, T) . (2) 

Here, the constant Wq stands to characterize the strength of scattering. We assume a 
glassy host and take the spectral density function in the Ohmic form S{Eji — E^f) = 
\Ej3 — Eal/tcG widely used in the theory of dissipative systems (see, e.g. Ref. [32]). The 
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temperature T enters into this expression through the function F{Ep — E^, T) defined 

as 

'^^'''-''-''^-[niE^-E,), E,<E^ ' 

where n{Ep — Ea) = [exp(\Ej3 — Ea\/T) — 1]^^ is the occupation number of the vibration 
mode with frequency \Ej3 — E^l/h. The term 

N 

^o.,^ E i:\<n?mj' (4) 
s=G,Cn=l 

represents the overlap integral of electronic probabilities for the states ipa and if) p. 

We describe the process of charge transport by means of the Pauli master equation 
for the populations Pa of the eigenstates a: 

2N 

Pa=ri{ft-Pa)+T^{f^-Pa) + J2 [{l-Pa)W^f,PfS-il-P^)WpaPa\ (5) 

(3=1 

where a = 1, 2, ■ ■ ■ , 2N and f^'^ are the Fermi distribution functions for the left and 
right leads: 

fa ' = [1 + exp jj , (6) 

fiL = Ep + eV and ^ir = Ep are the chemical potentials of the left and right leads, Ep 
is the Fermi energy at equilibrium taken to be in the middle of the non-disordered DNA 
band gap, which is the case for Au contacts [33]. The terms = 7(|'?/'f P + I'?/']" and 
= 7(IV'^P + IV'^P) measure the coupling between leads and the eigenstate a, with 
the parameter 7 being the strength of the coupling. 

We are interested in the steady state solution of Eq. ([5]). Solving the corresponding 
system of nonlinear algebraic equations by an iterative method that guarantees the 
condition < Pq < 1, we obtain the stationary current as: 

2N 

m = E Ff (/i^ - Pa) . (7) 

a=l 



4. Results 



In all calculations, we considered 30 base-pair poly(G)-poly(C) DNA molecules and 
used the following model parameters: unperturbed site energies = 1.14 eV and 
= — 1.06eV [34], while hopping integrals were adjusted [22] to reproduce the current- 
voltage characteristics measured in experiments on dry poly(G)-poly(C) DNA [5]: 
tcG = tec = 0.27 eV and tec = 0.25 eV. These values are within reasonable parameter 
intervals [35]. We assumed that the temperature is slightly below the freezing point of 
the environment (T = 273K) which allowed us to neglect all effect related to dynamic 
disorder (time-dependent fluctuations of the conflguration of solute ions surrounding the 
molecule) and validates the model, namely we consider only interaction with the phonons 
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Figure 2. Current- voltage characteristics obtained for different magnitudes of disorder 
indicated in the legend (upper panel) and the corresponding density of states (lower 
panel). 

of the thermal bath. The parameter 7 was found to only influence the amplitude of the 
current and was taken to be 7 = Wq in all calculations. 

In the upper panel of Fig. [2]we present the calculated current- voltage characteristics 
obtained for different magnitudes of disorder A = 1.5 eV (indicated in the legend). 
The lower panel shows the corresponding density of states. These results were averaged 
over 100 realizations of disorder for each value of A. The figure demonstrates that the 
electronic gap in the density of states is shrinking upon increasing the disorder while 
the voltage gap becomes wider, which could seem counterintuitive. The closure of the 
electronic gap is the usual consequence of disorder. It is due to the appearance of 
disorder-induced states in the bare band gap, which form the tails of the density of 
states. It may seem that these states can give rise to current at smaller voltage drops 
and to the consequent shrinking of the voltage gap. Fig. [2] shows just the opposite 
dependence of the voltage gap, namely it broadens with increasing disorder. 

In order to understand this counterintuitive dependence of the voltage gap, the 
following reasoning is in place. Although there are states in the energy gap (i.e. between 
bare band edges: \E\ < 0.5 eV in Fig. [2]), these states are strongly localized, so that 
direct tunneling via them is suppressed. In this case, only phonon-assisted hopping 
can give rise to charge transport. If two states are localized at different segments of 
the DNA, their wavefunctions have exponentially small overlap X^^ [see Eq. (jl])], and 
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2477 ' 



(8) 



therefore, the scattering rate W^p between the two states, given by Eq. ([2]), is small. On 
the other hand, for two well-overlapped states, the typical energy spacing is large due 
to the quantum mechanical level repulsion [36]. Such level spacing grows on increasing 
the disorder strength, reducing the thermally activated conductivity. 

The typical energy separation 5E between well-overlapped states can be estimated 
using the argumentation described in Refs. [37-39], where it was applied to the energy 
range close to the band-edge. Here, we use a similar reasoning for the band center. 
Briefly, a state extended over a segment of size N* senses not the bare on-site disorder, 
but rather a reduced one (averaged over the localization length) with a typical magnitude 
of a /\/ N*, were a is the standard deviation of the on-site disorder distribution (A/vT2 
in our case). The energy separation 6E is a function of A^* (the latter magnitude is 
considered to be large: A^* > 1). Then the equality 5E{N*) = A/\/l2N* provides a 
self-consistent estimate for 6E{N*) (see Refs. [37-39] for more details). Applying this 
equality to states in the center of the band associated with the G strand and taking into 
account that, for these states, SE{N*) ^ 27r|tGG|/-^* we obtain 

A 

^GG 

For the typical considered disorder, A 1 eV the above estimate gives 6E ^ 50 meV, 
the energy being larger than the considered temperature (~ 25 meV). The latter means 
that the phonon occupation number in ([3]) is small: n{\Ep — E^l/T) <^ 1, and suggests 
that phonon-assisted scattering from lower to higher states is suppressed and almost 
no temperature activated transport can take place in the system; hops from higher 
to lower well- overlapped states constitute therefore the dominant scattering process. 
Then, the only way for a charge carrier to hop from one lead to another is to make a 
series of cascade-like down hops over well-overlapped states. Appropriate sets of such 
states can only be found in the energy region above the bare energy band edge [31,40]. 
Disorder smears out the band edge, pushing the boundary of this region of "conducting" 
states up towards the band center (this boundary is analogous to the diffusion mobility 
edge [40]). Thus, for a larger value of disorder, such diffusion mobility edge lies at a 
higher energy and, therefore, a greater voltage is required to induce electric current 
through the system. The latter explains the observed dependence of the voltage gap on 
the disorder magnitude. 

The number of cascade states is of the order of N/N*, where N* stands for the 
localization length in the appropriate part of the spectrum. Because consecutive cascade 
states should overlap well, they are separated by the energy of the order of 6E given by 
Eq. ([8]). Then the smallest difference between the highest and the lowest energy levels 
in a cascade-like set of states can be estimated as 

(9) 



N iV , , 

AEjnin ~ ^ ; ^GG 



^GG 



In deriving of Eq. (Q we have used 6E ^ A/vT2iV* (see the preceding paragraph for 
details). Bearing in mind that the lowest state in a cascade should be about the bare 
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band edge Eq ^ 0.5 eV, we conclude that the quantity -Eo + '^-E'min shows the dependence 
of the voltage gap on both the disorder and the system size. The former supports the 
above argument on the increase of the gap with disorder while the latter suggests that 
the voltage gap increases with the system size. 

Note that the latter trend could be reversed for very short DNA chains whose length 
is of the order of the smallest typical localization length. In this situation all states are 
coupled to both leads and contribute to charge transport. Thus, the current appear 
as soon as the Fermi level of a contact aligns with the lowest state in the energy gap 
and therefore the voltage gap coincides with the energy gap and decreases with disorder 
(see Fig. [2]). However, in experiments the chain length is typically much larger that the 
charge carrier localization length, ^ A^*, so the dependence given by Eq. ([9]) should 
hold, which suggests, in particular, that long DNA chains are insulating. 

5. Conclusions 

We considered theoretically the charge transport through a synthetic double-stranded 
poly(G)-poly(C) DNA molecule attached to two ideal leads and embedded into a random 
static environment (e.g. solvent below its freezing temperature). A nonlinear Pauli 
master equation for the populations of localized electronic states was used to describe 
the hopping transport of charge carriers. We demonstrated that the voltage gap becomes 
wider as disorder increases. The calculated I — V curves could indicate that the 
conductivity is band-like, although the charge transport is incoherent. This suggests 
that care should be taken when concluding on the nature of the charge transport 
by examining only I — V characteristics. The proposed method of electric current 
calculation is applicable for a broader range of systems, e.g., organic polymers. 
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